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We study the Casimir effect due to a massive vector field in a system of two parallel plates made of 

real materials, in an arbitrary magnetodielectric background. The plane waves satisfying the Proca 

equations are classified into transverse modes and longitudinal modes which have different dispersion 

relations. Transverse modes are further divided into type I and type II corresponding to TE and 

TM modes in the massless case. For general magnetodielectric media, we argue that the correct 

boundary conditions are the continuities of Hy , (j), A and dxA^, where x is the direction normal to 

the plates. Whereas there are type I transverse modes that satisfy all the boundary conditions, it 

is impossible to find type II transverse modes or longitudinal modes that satisfy all the boundary 

f^ ' conditions. To circumvent this problem, type II transverse modes and longitudinal modes have to 

fSJ ' be considered together. We call the contribution to the Casimir energy from type I transverse modes 

_. , as TE contribution, and the contribution from the superposition of type II transverse modes and 

ijT^' longitudinal modes as TM contribution. Their massless limits give respectively the TE and TM 



contributions to the Casimir energy of a massless vector field. The limit where the plates become 
perfectly conducting is discussed in detail. For the special case where the background has unity 
C"*^ ■ refractive index, it is shown that the TM contribution to the Casimir energy can be written as a 

sum of contributions from two different types of modes, corresponding to type 2 discrete modes 
and type 3 continuum modes discussed by Barton and Dombey ijl§|. For general background, this 
splitting does not work. The limit where both plates become infinitely permeable and the limit 
where one plate becomes perfectly conducting and one plate becomes infinitely permeable are also 
investigated. 
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C^ ■ I. INTRODUCTION 

en 



In 1948, Casimir [1| predicted the existence of an attractive force between two perfectly conducting plates as a 

manifestation of the zero-point vacuum fluctuations of quantum fields. Although experimental veriflcation of this 

prediction has little progress before 1997, Casimir effect has attracted more and more attention since 1970s 2]. It was 

gradually realized that this effect is closely related to other areas of physics such as quantum held theory, cosmology 

and condensed matter physics. More recently, the advances in the Casimir measurements and nanotechnology have 

• • . stimulated more interest in Casimir effect. 

.5^ ' Historically, the Casimir effect was investigated for electromagnetic field (massless vector field). In 1956, Lifshitz 

^^ , [3| derived a formula to express the Casimir force acting between two dielectric semi-infinite slabs in terms of the 

;_( ' permittivity of the slabs [^Tal • Refinements and generalizations of the Lifshitz formula where the semi- infinite dielectric 

5r '• slabs are replaced by dielectric or magnetodielectric plates of finite thicknesses were derived in [9l-[lll [iSl - fm [23j . In the 

limits the permittivity goes to infinity, Lifshitz formula reproduces the result of Casimir for two perfectly conducting 

plates. 

Since 1970s, Casimir effect was studied for other quantum fields such as scalar fields and spinor fields, for both 
massless and massive sectors. In contrast, there are very few works that studied Casimir effect of massive vector 
fields. One of the incentives to study Casimir effect of massive fields comes from the prevalence of spacetime with 
extra dimensions which were proposed to solve some fundamental problems in physics. To study a quantum field in 
a spacetime with extra dimensions, one can use Kaluza-Klein decomposition to decompose the quantum field to an 
infinite tower of massive fields in AD. This approach has been intrinsically used in a number of works to study the 
Casimir effect of scalar field or spinor field in spacetime with extra dimensions such as Kaluza-Klein spacetime and 
Randall-Sundrum spacetime. For electromagnetic field, the Casimir effect on a pair of perfectly conducting plates in 
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the Kaluza-Klein spacetime M^ x S^ was studied in [17] using the results of Barton and Dombey [T^ on Casimir 
effect of massive photons. 

The pioneering work of Barton and Dombey [18l.[l9| considered the Casimir effect on a pair of perfectly conducting 
plates due to a massive vector field. Unlike the scalar and spinor case, the extension from massless vector field to 
massive vector field is in fact much more complicated. The classical Maxwell's equations have to be replaced by Proca 
equations |20j and the gauge degree of freedom is lost. Therefore, a massless vector field has two polarizations but 
a massive vector field has three. For a pair of infinitely large perfectly conducting plates, the contribution to the 
Casimir energy from the two polarizations are the same. However, the contributions to the Casimir energy from the 
three polarizations are not the same. Another complication in generalizing Casimir effect of massless vector field to 
massive vector field is that besides the conventional boundary conditions which require the continuities of Dj^, Eii , Bj^ 
and Hii across interfaces, one also need to impose the continuities of the scalar potential (j) and the vector potential 
A [21|. For massless vector field where we usually impose the gauge conditions ^ = and V • A = 0, A_l is not 
continuous for the farnilv of TM modes. This gives an obstacle to generalizing the results from the massless case to 
the massive case. In [1^, Barton and Dombey showed that the eigenmodes of a massive vector field can be divided 
into two families of discrete modes and one family of continuum modes. In the massless limit, the contribution to the 
Casimir energy from the discrete modes yields the result of Casimir [Ij- The contribution to the Casimir energy from 
the continuum modes has zero massless limit. One thing that worth remarked is that one family of the discrete modes 
corresponds to the TE modes in the massless case, but the other family of discrete modes is not a direct generalization 
of the TM modes in the massless case. 

As a first step to investigate the Casimir effect on a pair of plates made of real materials in a spacetime with extra 
dimensions, we generalize the work of Barton and Dombey [T^ to a pair of plates made of real materials in some 
magnetodielectric background medium. This can also be considered as generalizing the Lifshitz formula for massless 
vector field to a formula for massive vector field. As explained above, the generalizing of the Casimir effect from 
massless vector field to massive vector field is nontrivial even for a pair of perfectly conducting plates. Therefore, 
one should expect that the result will even be more complicated when one considers real materials. As was discussed 
in 18], for dielectric plates, the plane waves of a massive vector field can be divided into two families of transverse 
waves and one family of longitudinal waves. However, one of the families of discrete modes for perfectly conducting 
plates is in fact not transverse, and the continuum modes for perfectly conducting plates is in fact not longitudinal. 
Therefore, it is not obvious how one should generalize the work of Barton and Dombey [18|. In fact, one can show 
that for plates made of real materials in a magnetodielectric background, only one family of the transverse modes 
can satisfy all the boundary conditions. As a result, to look for other eigenmodes of the massive vector field, we are 
forced to consider the superposition of the other family of transverse modes and the family of longitudinal modes. 

In [3|, Lifshitz modeled two dielectric plates as two semi- infinite slabs which is a reasonable model when the skin 
depths of the materials of the plates are much smaller than the thicknesses of the plates. The configuration of two 
semi-infinite slabs separated by a medium is a three-layer model. However, since we want to consider plates made of 
any materials, we find it necessary not to make any assumption about the skin depths. Therefore, we would consider 
two plates with finite thicknesses which can be modeled by a five-layer configuration. 

The layout of this article is as follows. In Section |lTl we review the Proca equations for massive vector fields. In 
Section IIIIl we divide the plane waves solutions to the Proca equations into type I and type II transverse waves, 
and longitudinal waves. In Section lIVl we discuss the boundary conditions that should be satisfied by the potentials 
and fields on the interface of two plane parallel media. In Section |Vl we derive the Casimir energies contributed by 
transverse modes of type I, and by the superposition of transverse modes of type II and longitudinal modes. We show 
that their massless limits are the TE and TM contributions to the Casimir energy due to a massless vector field. In 
Section IVTl we discuss the limiting case where the plates become perfectly conducting. In Section rVIIi we consider 
the limit when the plates become infinitely permeable. In Section [Villi we study the case where one plate is perfectly 
conducting and one plate is infinitely permeable. 

II. PROCA EQUATIONS OF MASSIVE VECTOR FIELD IN A MAGNETODIELECTRIC MEDIUM 

A massive vector field is represented by a four- vector ( — , A^ , Ay , Az), where 4> is the scalar potential and A = 

{Ax,Ay,Az) is the vector potential. Define the electric field E and the magnetic field B by 

a A 
E = — ^-V0, B = VxA. (1) 

dt ^ ' 

Consider a magnetodielectric medium with permittivity e and permeability /i. Let pf and 3 f be the free charge and 
free current in the medium. Assume the usual linear relations D = eE and B = /iH. Then the Proca equations for 



the electromagnetic waves propagating in this medium are |18| : 

V • B = 0, (2) 

VxE+ — =0, (3) 

9 

The first two are the well-known Maxwell's equations which are automatically satisfied because of ([TJ. The conserva- 
tion of free charges implies the continuity equation: 

Differentiate dH) with respect to t and apply the divergence operator to ([5]), (O implies that the Lorentz condition 



f 
? dt 



:^^ + v.A^o (7) 



has to be satisfied. From ([T]) and (O, one can then derive from (|4|) and ^ the following two equations for (p and A: 

^M^^ - ^' + -^ j A - {e^l^ - 1) V(V • A) =mJ/. (9) 

([7]), ([5]) and © are the equivalences of the Proca equations for the potentials (j) and A. 

For a massless vector field, the primary quantities in the Maxwell's equations are the electric and magnetic fields. 
There is a gauge degree of freedom given by 

A^ K + Vip, (h^(j)-^ 

at 

for an arbitrary function ip. The Lorentz gauge ([7]) is one of the gauge conditions that can be used to fix the gauge. 
This gauge degree of freedom is lost in the massive case. In fact, for massive vector field, the primary quantities are 
the potentials 4> and A, and the electric and magnetic fields are derived quantities. The Lorentz condition ([7]) is a 
necessary condition followed from the conservation of charges. 

III. PLANE WAVES AND DISPERSIONLESS RELATIONS 

Consider an unbounded magnetodielectric medium with permittivity e and permeability p. in the absence of free 
chargers and currents, i.e., p/ = and J/ = 0. As discussed in [18| . the monochromatic plane wave solutions of the 
Proca equations can be divided into transverse waves with V • A = and longitudinal waves with V x A = 0. For 
the transverse waves, it follows from Lorentz condition (O that = 0. The equation ([9]) for the vector potential A 
becomes 

e,^-V^ + ^JA = 0. (10) 

Recall that for a massless vector field, the vector potential A satisfies the equation P^ with m = 0, and the conditions 
(j) = Q and V • A = 0, called Coulomb gauge or radiation gauge, are usually imposed to remove the gauge degree of 
freedom. Therefore the transverse waves are in one-to-one correspondence with the massless electromagnetic waves. 
Since massless electromagnetic waves are usually divided into TE and TM polarizations, we do the same for the 
massive case, but call them type I and type I I transv erse waves. 



In the following, let k^ — (/c2,A;3), k± = \/k\ + fc| and define f'k±,u:{y, z,t) 
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A. Type I transverse waves. The type I transverse waves are the equivalence of the TE waves in the massless 
case: 






he'p^-h^Ay,^,t) 



In this case, E^ — 0. The dispersion relation is: 

2 2 
2 2 J 2 ^^ ^ 



tf 



B. Type II transverse waves. The type II transverse waves are the equivalence of the TM waves in the massless 
case: 

Pt 
In this case, B^ ~ 0. The dispersion relation is: 

2 2 



fi2 



C. Longitudinal waves. The longitudinal waves are 






The dispersion relation is: 



2 2 



= 0. 



For longitudinal waves, V x A = implies that the magnetic field vanishes identically, i.e., B = 0. On the other 
hand, the electric field is given by 



E^, 



E,, 



E. 
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pi-fcl e'P-^/k,,^(y,z,t) 



ifcsc / UJ 



to 



— -pi-fci e^P-Vk,,.(y,z,t) 



tK2m c^ 
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e'p-^h^Ay.^.t) 



e'p^-h^Ay,z,t) 



Notice that the electric field also vanishes in the massless vn — >■ limit. It should be observed that the longitudinal 
waves and transverse waves satisfy different dispersion relations, unless e/i = l/c^, i.e., the waves traveled at the speed 
of light in the medium. 



IV. BOUNDARY CONDITIONS 



In this section, we consider the boundary conditions that must be satisfied by the potentials A and (j) and the 
fields E and B on the interface of two plane parallel media in the absence of free charges and currents. Let x be the 
direction normal to the boundary of the two media. As in the massless case, © and ^ imply that B^ = B^ and 
E|| = {Ey.E^) are continuous. Taking the derivative of (Hj) with respect to t gives 



(dTi\ m^c^ I d(j) _ 

' \dt) ^"^c2 W~ ■ 



The Lorentz condition ([7]) then imply that 






A =: „. T^At 



Therefore, we find that 

('- 

\ dt ^ih? ) ^ dt fih? 

must be continuous. ([5]) then implies that Hy = — -, — - ] has to be continuous. Notice that the boundary condition 
that one obtains from ^ is the continuity of 

dT> m'^c^ . 
A 



dt nft^ 

but not the contin uity of D^ as in the massless case. 

As discussed in [18J, the boundary conditions above are not enough for massive vector fields. It was pointed out by 
Kroll [21] that the potentials <j) and A also have to be continuous. Therefore, the complete set of boundary conditions 

is the continuities of Bx:Ey,Ez, — tt— ^^x, — -, — -,Ax,Ay,Az,(j). However, not all these conditions are 

dt ^h^ M A^ 

independent. For example, the continuities of and Aii — {Ay,Az) imply the continuity of 

The continuities of — ^, — - imply the continuity of — —-^ pr^a; since ([S]) gives 

[i, [i. at [ihr 

d{eE^) w?^ 1 fdBz dB. 



Ax — — 



■"y 



dt fih? fi \ dy dz 

Similarly, ([3]) implies that if E|| = {Ey,Ez) is continuous, then so is i?^,. On the other hand, the continuities of 
A|| — {Ay,Az) and the Lorentz condition 

1 50 , 1 30 OA.^ 

dA^ . 
imply that (j) is continuous if and only if -- — is continuous. 

dx 

In conclusion, for the boundary conditions, it is sufficient to impose the continuities of A^, Ay, Az, — -, — - and the 

/i ^^ 

continuity of cither (j) or —- — . 

ox 

V. CASIMIR ENERGY OF PARALLEL MAGNETODIELECTRIC PLATES INSIDE 

MAGNETODIELECTRIC MEDIUM 

In this article, we are interested in computing the Casimir energy of two parallel magnetodielectric plates inside a 
magnetodielectric medium (see FIG.[T]). We assume that the cross section of the plates are infinite. For convenience. 
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FIG. 1: Two parallel magnetodielectric plates inside a magnetodielectric medium. 




FIG. 2: A five-layer model 



we first consider a five-layer model consists of five plane-parallel layers of magnetodielectric media as shown in FIG. 
[21 At the end, we only consider the case where ei — e^ — £5 and ^.i — jjlj, — ^q. The interfaces of the media are 
located at x = ai,a2,a3 and 04, with ai < a2 < 02 < 04. We assume that there are two artificial boundaries at 
X = flo < ai and x = 05 > a4, and let di = Ui — a,;_i, 1 < « < 5. At the end, we are going to let di, ^5 — >■ 00. We also 
set ti = d2 and tr = d^ to be the thicknesses of the left and right plates that we are interested in, and set a = 1^3 to 
be the separation between the plates. 

For the Casimir energy of massless electromagnetic field in the five layer model, one can consider the contribution 
from TE modes and TM modes separately. However, it turns out that for the massive case, one cannot consider the 
contributions from the three types of waves separately. In fact, the contribution from transverse modes of type I can 
be considered independently, but the contributions from the transverse modes of type II and the longitudinal modes 
have to be considered together. It is impossible to find type II transverse modes or longitudinal modes that satisfy all 
the boundary conditions. 



A. Contribution from transverse modes of type I (TE modes) 



For transverse modes of type I, assume that 



A7, 



fc3(2t,e^P"-'^ + 'B,e-'P--^)/k,,^(2;,z,i) 



(h=0 



for ai_i < X < Oi. 



The dispersion relation implies that 



The electric field E and the magnetic field B are given respectively by 

E, -iL.fc2(a,e'P--^ + <8.e-'P--^)/k,,^(y, z,t) [ B, = ^pT,^h{Qi,e'P^-- - »,e-*''--^)/k,.^(y, z, t) 

for Oi-i < a; < a^. In this case, 4>, A^ are automatically continuous. The continuities of Ay,Az imply that for 
1 = 1,2,3,4, 

The continuities of — -, — - imply that for i = 1, 2, 3, 4, 

Next, we have to impose some boundary conditions on the artificial boundaries x = a^ and x = a^. As in the massless 
case, we impose the conditions Ey = Ez = and B^ = 0, which gives 

From here, one can use the method we developed in our previous work 22]. After a renormalization by the criterion 
that the Casimir energy should vanish in the limit the plate separation ^3 is infinite, and passing to the limits 
di, (is — ?> 00, one finds that the TE contribution to the Casimir energy is given by 



■ r /'°°ln|l-D™(^,fc^)e-25T,.(?.fc^)aU^^fc^^^^ 



Ec^s~ I I ln<!l-D'^'-(e,fc^)e-^«--(^--^'^^fc^dfc^de, (11) 



where 



A^Al (^ — p-lT,lti\ /I _ p-qT.rtr\ 

s™(e, fcx) = , ""[^^X x/^^ . — ^- (12) 



(1 - [A|]^e-9^.'*') (1 - [Ai]^e-«^>-*-) 



Here the variables have been renamed by ei = £3 = £5 = £;,, /ii = /is = /is = /ib, £; = £2, Sr = ^4, Mi = M2, Mr = M4> 
ti = d2,tr = d4,a = ds] and 



(13) 



The TE contribution to the Casimir force is given by 

(5£:TE ft /-oo /.oo 






:7^/ / <ZT,6(e,fc±)([S™(e,fc±)]''e'^^'^(«''-^'^-l) k^dk^dt 

^TT Jo Jo ^ ^ 



Obviously, in the massless rn — )■ limit, one obtains the result for the massless case, where the mass m in qxji^, k±) 
in (fT51) is set to zero. 



B. Contributions from superposition of transverse modes of type II and longitudinal modes (TM modes) 



For the superposition of transverse modes of type II and longitudinal modes, we have for a^-i < x < ai 



■^x — 

A, 



li(e:^ew.x _ S),e-'P^-*") +pL,»((£,e*P^-" - ^^e-'P^-'') 



PTa 



h^.Ay^z,t) 






(14) 



The electric field E and the magnetic field B are given respectively by 



E^: =i 



Ey =i 



E,^i 



jk 



T)t Tfl C 



Ei^iiuh? 



PT,, 

J 2 2 

JcoTn c 



By ^^ (p|,, + fci) (£,e^P--- - S),e-^P---)/k,,^(y,z,t) 
S. = - — (pt . + A:i) (C^e^P--^ - S).e-*P---)/k,,„(y, z, t) 

PT,i 



h^,uj{y,z,t) 
h^,uj{y,z,t) 



(15) 



(16) 



Since B^ = 0, we can call these modes TM modes. The continuities of — -, — -, (p, A^, Ay,Az imply that for i — 1,2, 3, 4, 



1 



{Ph + kl) (C^e^P--"' - D,e-^P--"*) = 



1 



{Ph+i + ^i) (2:^+16^^--+^'^' - 2).+ie 



-ipr.i+ia; 



fJ-iPT,i ' lH+lPT,i+l 






,iPL,iai 



die 



-iPL,iai 



(17) 



-(e:,+ie^P^.*+i"' - S,+ie-^P^-*+i'^') +PL,^+l(*S^+le 



,«Pi,i+iai 



i?j+ie' 



-ipL,i+iai 



PT,i+l 

- (e:,+ie^P^.*+i"' +S),+ie-'P^-+i"0 - (e»+ie*P^-+i'^' +^,+ie-^P^.'+i'^')- 

For the artificial boundaries at x = ao and x — a^, one can impose the conditions Ay = A^ = (j) = 0. These imply 
that 



£56*^^'^°= + 2)5e-'P^-="s =0, ^56*^^'="= + ^ge^'P^'^"^ = 0. 
As in the case of TE modes, one can then show that the TM contribution to the Casimir energy is given by 



pTM _ 



^0 



In det — — 7-^ — ; — -dS^k^dkj 

Qoo(4,fc-L) 



(18) 



With the notations 

qL,j{S.,k±) 



fc2 



SjiiO^^iiiOfi'^ 



, J "^^b, I, r; 



_ A^fc(^OgT,b(C,fc±)[gT,j(^,fc±)^ - fej] - Mj(^C)gT,j(^,fc±)[gT,fc(^,fci,)^ - fcj 

Aifc(«OgT,6(C,fc±)[gTj(^,fc±)2 - A:^^] +Mj(*C)gT,j(^,fc±)[gT,fc(^,fc±)2 _ ^^2^ 



■, i^l,r; 



^1 l?i'^±J— T^ TTT-;;^ — V77 TT^T^ — v5 7_2 i , — 7 TT^: — T7I TT^T — v5 TTT' J— 'I'j 



9Lj(e,fc±)bL,b(e,fc±)2 - fci] + gL,;,(e,fc±)[<ZL,j(e,ft±)' - fci 



"jU2(C,fc±) = 



9LJi(^,fc±)^ -qL,J2itk±f 



'iL.n{tki_Y-kl 



ji,J2 =b,l,r; 



ftlj2(C,fc±) =- 



9Tji(C,fc±)'ZLj2(C,fc±) 



. Mj2(»0 [gT,j-i(C,fc±) -fc±] A . . 



= fc,/, 



+ f^ fc , ) =1 + gr.J2 (C, k±)fij, (iQ [gT,ji (g, fc±) ~ fc^] _ 



'*iU2'-^''*^i 



is a 4 X 4 matrix with components given by 



, ji,J2^b,l,r; 



^1 = 



— A}^ (l — e -1^.'-' ) e'' 



g-2qL,lti\ ^qL,lti 



'lb ' bl 






^13 






-qT,ba 



= (l-[A:,i]%-2,T,.t.Ag7T,.t. 






.14 = — — (1 - Ajie-29--*-) e-?--*"- + ^ (1 + Anie-29.,.*.) g?...*.^ 



K^ 



Q21 = (1 - [Aji] ' e-29T,*A ^,T,h + ^^ (1 _ 



,-2qi^jti\ qL,ltl 



Q2 



^bi ^ib 



aih 



(1 - A|ie-29T,,t,) gg^,,t, + ^ (1 + A|"e-2,^.'*') e«^-'*S 



g23 = 



/24 — 



:?31 



/32 — 



-A^^ (1 - ^ 



-2qT,rtr\ pQT.rtr i ^^bRbr 



-^ (1 _ e-29.,.*.) g«.. 



g-qT,ba 



_^!xb (All _ g-2g,,.t.) ^qr^^u + ^ (1 + A^e-^?^-*-) e"--*^ 

'^hr '^rb ^rb 

_Plb /-j^ _ ^ll^-2qT,iti\ ^qr.lh _ ^^ ^^III _|_ g-29i,,t,\ g<?i,it, 
'lb ' bl ' lb 

ablPlb 

Hl'^lb 



-qL.bO- 



(1 - e^^-jT.it^^ g<?T,,t! _ A|" (1 - e"2<?i,jti) g9L 



)33 =^ (1 _ Aiie-29T,.t.) g9T,.t. + ^^ (1 + A™.-2,,,.t.^ ^9,,.*. 



' rb 
PrbCtbr 



Pbr Kr 



br rb 



g34 ^Bip^ (1 _ g-2«T,.*.) g?T,.t. + ^ [Aj,"]'e-29.,.t.A g,.,.t. 

^rb^br 



i,!*!"^ plL.ltl 



' bl 'lb 
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g42 =^1^ (1 _ g-2,T,*,) ggT,t, ^ ^ _ [^III] - g-2,,,*, gg,, 



Hl'^lb 



/43 — 



/44 



-^ (1 _ ^ll^-2qT.rU^ ^qT..U „ ^^b (^III ^ g-2g,,.t.^ gg,,.t. 
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and 



/ Qi3 Qi4^ 

Q21 Q22 

Q33 Q34 

VQ41 Q42 

Q . 
is obtained by taking the limit a — > 00 of Q. In general, det -- — is a very complicated function of the distance a 

^; 00 

between the plates. For the Casimir force acting between the plates, let Mij be the minor matrix of Q obtained by 
deleting the i^'^-row and j^'^-column from Q. Then the TM contribution to the Casimir force acting between the plates 
is given by 

TM_ dE™ _ h [^ f°° f qrAQiiMii - Q12M12 - Q23M23 + Q24M24) 



CiilVl _ _ _ 

^''' da 47r2 J^ J„ { det Q 

, gL.b (Q31M3I - Q32M32 - Q43M43 + Q44M44) 1 ,^ ,, ,. 

+ d^tQ I fc±dfc±d^. 

In the massless m ^ limit, qi^b = QL.i = QL.r- Therefore, aj-^j^ — for all ji,J2 — b,l,r, and Aj^^ — A^}^ — 0. 
The matrices Q and Qoo reduce to matrices of the form 




Ql- 







\Q4l Q42 




Consequently, 



1 , Q _ _ Q42Q34 \Q11Q23 Q21Q13) _ -, _ Q11Q23 

"'^ 7to~ ^ n" n" n" n" ^ 7W7W 

^/oc V42V34V21V13 



34^21^13 ^21^13 

A" A" (1 - e-29T,,tA (1 _ g-2,^,.t.) (19) 

where the massless limits of Aj\j — l,r, are 

hm A" = A^fcgT,b(£jMj$^) - MjgT,j-(£bA^fc^^) ^ gT,6(g,fc_L)L^o£j(ig) - gTj(C, fc±)L^o gb(^C) ,20) 

One finds that the massless limit of the contribution to the Casimir energy from the superposition of the transverse 
modes of type II and longitudinal modes is precisely the contribution from the TM modes in the massless case. This 
is another justification that we call the contribution from the superposition of the transverse modes of type II and 
longitudinal modes as TM contribution. It is interesting to note that in the massless limit, the two polarizations for 
the massive field naturally reduce to only one polarization for the massless field. 

VI. THE CASIMIR EFFECT ON A PAIR OF PERFECTLY CONDUCTING PLATES 

In this section, we study the limits of the Casimir energy and Casimir force when the plates become perfectly 
conducting. It is well-known that this can be obtained by the limit £; = e^ — > 00. In this limit, we find that 



qTj,qT,r — > 00; qr^h = V£fc(«C)Mb(«C)C^ + k 



h? 




IL., -M- + ^^ + sSOl^SiW '"'' ^ '"" -^ V - + /^^ - 'Zo, 
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90 gi,fc - fci + 9l,6 [ql - k'i) 



'^br'^^br 



k\ 




9o - 9i,6 




qr.b 


qh.b{qh - 


-fcD + ^olgi^fc- 


-fci)" 



OilbPlb OLuPbl OLrbPrb CUbrPbr , ■-^ :.>j ^n,u . . 



For the TE contribution, it is easy to see from ([T^ that T)'^^ -^ 1, and thus the TE contribution to the Casimir 
energy (fTTj) reduces to 



^™ -^ [ [ In (l - e-^'-^(«-'=-)'^) fc^dfcxd^ (21) 

The TE contribution to the Casimir force is then given by 

F^i =-^,r r '^^m, M dk.dt (22) 

^T^ Jo Jo exp (25T,b(C, k_L)a) - 1 

which is always attractive. It is interesting to note that in the perfect conductor limit, the TE contribution to the 
Casimir force is independent of the thicknesses of the plates. 

The TM contribution is much more complicated. With the help of computer, one can show that in the perfect 
conductor hmit, the TM contribution to the Casimir energy (1181) is given by 



where 



pTM 6 r r.. rni.kj 



"---ii^l I '-ndpn'-*-'^' <===" 



W^ ^ ([1 - A]2 - [A + A]2e-2*t,) ([1 _ ^]2 _ [A + AJ^e-^wt,.) ^ 



and 



W =Woo - ([1 - A2] - [A2 - A2]e-29ot,) ([1 _ a2] _ [a2 _ A2]e — ) e 



([1 - A] [A - A] - [1 + A] [A + A]e-29oti) ([i - A] [A - A] - [1 + A] [A + Ale^^-Jot.) g-2gi„,a 



(24) 
+ 4A(1 - A) ([1 - A] + [A + A]e-29ot,) ([i _ A] + [A + AJe-^-Jot.) g-?T,6ag~gi„,a 

+ (1 - A)2(A + A)2 (l - e"290t,^ (-j^ _ g-2got.^ ^-2qT.,a^-2qr.,ta^ 

For an arbitrary background medium, one cannot spht the Casimir energy psp into contributions from two different 
polarizations. In the special case when the background material has refractive index nt = c-JShUb equal to one, e.g., 
if the perfectly conducting plates is placed in vacuum, then 



C^ , 7 2 I "^^'^^ 



qT.b = qL,b == V ^ + ^i + -^T' •= '2'™' 



K^ 



and W/Woo simplifies drastically to 



W ^(. g-2,„.^ (, (1 - A)2(A + A)2 (1 - e-2.0^-) (l - e-^^Q*-) _^^J 

\ ' \ (Tl _ A12 _ FA ^ ^^2^-2Qnt,\ m _ AI2 _ fa ^ A12^-2ant^A 



M^oo ^ M ([l-A]2-[A + A]2e-29ot,)([i_A]2-[A + A]2e-29ot.) 

/ jj2(i_g-2,o.)(i,g-2,o..) ^^^^^^ 



(l_e-29„a) 1 



where 



(1 - D2e-2qotl) (1 _ £)2g-2got,.) 



^ A + A 



1-A 

,TM, I J T^TM, II 



In this case, the TM contribution to the Casimir energy (|23p can be written as a sum of two terms £^ca,s' ^^^^ -^c 

^Cas 



-^Cas' i^ ^^^ same as the TE contribution, obtained from (l?T|) by setting ri;, = 1, and it corresponds to the massless 
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TM contribution. The corresponding Casimir force is given by ((22)) with rib = 1, which is always attractive. The 
contribution to the Casimir energy from the second TM polarization is given by 



E, 



TM, II 
Cas 



h 



OO /"OO 



InO- 



D{^, k^f (1 - e-2«o*') (l - e-29«*-) 
(1 - D{i, k^ye-^iot' ) (1 - D{^, fci)2e-29ot.) ' 



-2«'"'^ I k±dk±d£,. 



(25) 



Some simple computations give 



Di^, fcj 



, 2 



qo 



1-A 



go {ql - ki) + qra (q'o ^1) ^ (^o' - ql)t *' 



90 



The contribution to the Casimir force from the second TM polarization is then given by 



F, 



TM, II 
Cas 



h 

2^ 



OO poo 



{1-D{^,k^)'e- 



■2?o*i 



(1-D{^,k^ye 



2„^2qoU 



„2qma 



Since < Z? < 1, 



1 - i?2p-2got 
Z? (1 - e-29oi) 



L»(C, fc±)2 (1 - e-29oti) (1 _ e-29ot.) j 

1 - L)2e-2got _£»(!_ e-2got^ (i _ /)) (i + i:)e-29ot) 



^(1 



-2qot 



D{l-e-^'}ot) 



k±dk±d^. 



>0. 



(26) 



-,TM, II 



Therefore the sign of Z^cas' i^ always negative, i.e., it is an attractive force. 

Comparing ([25]) to pT|) . we find that the second TM contribution (|25p can be identified as the TE contribution to 
the Casimir energy of a pair of dielectric plates due to a massless electromagnetic field, where the permittivity of the 
dielectric plates is 

^2^4 

e{u;) = 1 - 



ft2a;2- 



This fact was observed in [17|. In particular, in the limiting case of perjfectly conducting plates in vacuum, our general 
formula yields the result that was obtained by Barton and Dombey [18]. 

In order to gain more insight how the two TM polarizations split out in this special case, as was described in 
[18| , let us return to the expression for the TM modes (J14p . In the case of perfectly conducting plates in vacuum, 
PT,i = Pla = p for i = 1, 3, 5 and pt,2tPt,a — >■ oo. Let p^^i — pq for i = 2,4. Notice that 



^^p^ + kl + 



ft2 



= pl + kl. 



The first boundary condition in (|T7)) then implies that C2 = '^2 = £4 = 2)4 = 0. Since £2, £4 —^ 00, ([T5|l and (fTB|) then 
imply that both the electric field and the magnetic field vanish in the plates. In the vacuum, i = 1, 3, 5, the potentials 
are given by 






fc2 

p 



P<B, 



— S)» 
P 



pdi 



h±.u{y,z,t) 



- k2 [{€^ - £,) e'P^ 

A, = -k3 [{<^^ - <^^) e'p- 



(D,;-5,)e-^^1/k,,c.(2/,^,i) 



0=_(p2 + A:i)(g,e^P-+5,e-'P-)/k,,„(y,z,i) 



In the plates, i — 2,4, the potentials are given by 

r A^ =po ((£.6*''°" 



Ay =fc2 (€.e*P°^ + ^,e-*P°") /k,.^(y, z, t) 

A, =fc3 {€,€'""'' +S,e-'P"^) h^Ay,z,t) 

=u; (€,e'Po- + y,e-'P«") h^^^V, z, t) 



(27) 



(28) 
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For i — 1,3,5, renaming the variables by 



fc2 

p 

fc2 
P 



or equivalently. 






p^ 



c 



p 



,^2 



^1 , ^i 



2 ^i 






p 



p 



pz 



-T>' 



P 



p2 + fc2 



'UiJ U'i 



p2 + A:2 -. 
fc2 



The formulas for the modes in the vacuum (1271) then becomes 



pz 



Ay =fc2 (C^e^P- + S)^e-*P-) /k,,.(y, z, i) 



S)' 



p2 + kj^ 
P 



pz 






for i — 1,3,5; and the boundary conditions (fT7|) becomes 



jpoa; 



(29) 



which follows from the continuities of cj) and Ax, Ay. Here / = 2,4. For I = 2, j = 1 and i = 1 or j = 3 and i = 2. For 
I = 4, j = 3 and z = 3 or j = 5 and i = 4. Substituting the third equation into the first equation, one obtains 



ple'e^P^'+y'e-^P'^O^PoCe/e 



ipoat 



+ die-'P'"''), 



which just reflects the continuity of 
the system 



dAx 
dx 



The boundary conditions thus can be decomposed into two systems, one is 



p(€'e'P'^'+r,e-'^"0=Po(e 



dAx 



;e'' 



S'/e" 



-tpai 



Po (€ie'P°''' - die 



Woo-i 



(30) 



that reflects the continuities of A^ and ^-^ , together with the conditions <B\ e*^"" +d\ e"*^"" = 0, C'r.e^P'''^ +5Ce"*^''= = 

ox 
at the artificial boundaries x = a^ and a; = as; and the other one is the system 



er'-e^P^'+D'-e^^P"' =£r 



S'ie" 






((£;-e*'^'+^;-e-*P'^") 



(31) 



that reflect the continuities of Ay, A^, together with the conditions £'16*^°° + £)ie-*P"° = 0, £56*^°^ + D^je"*''"^ = 
at the artificial boundaries x = a^ and x = 05. In fact, pip determines (£',S)'),j = 1,3,5, from (2;;,S';)'^ — 2,4, but 
«i (ioes not determine uniquely. One has the freedom to add any solutions satisfying 

ir^.gipa, ^ ^^^.e-^P"^- ^0, j^ 1, 3, 5, 

corresponding to the trivial solution of (€',5^'),j — 1,3,5, and {<S.i,di),l = 2,4. The nontrivial solutions of 
[^j5^j)i.? = 1)3,5, and {€1,^1), I — 2,4, satisfying (POJ) correspond to the type 3 continuum modes discussed in 
[18| . For j = 1, 3, 5, choosing 

^J T.2 J' J Ti2"J' 

Po Po 
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FIG. 3: The dependence of the Casimir forces on the mass m when the background medium has refractive index 1 and 2. 



eq. (|3ip is satisfied, and we can fix the type 3 continuum modes to be 



A -^ 



{^^''^ + ^'^ 



' h^,uiiy,z,t) 



A, 



{€ie'P'^+S',e-^P^)h^.Jy,z,t) 



'pI 

Po 



where i = 1,3, 5 for a^-i < x < a^. For ai < a; < 02 and 03 < a; < 04, the modes are given by (j28p . The boundary 
conditions satisfied by these modes are given by ((30|) . It is easy to verify that the contribution of these modes to 
the Casimir energy is given by ideas' ■ ^O'" ^^^ trivial solution £^ = 3"^ = 0, j = 1,3,5; £; = ^i — Q,l = 2,4, the 
potentials vanish identically on the plates. In the vacuum, we have three systems of equations: 



g-/gipaj_i 



j,/g-«pa,„i 







e^'.g«P<ij ^ J)'.g-*Paj 







1,3,5. 



These three systems are independent. In the limit di = ai — oq and ds = as — 04 go to infinity, we only need to be 
concerned with the system with j = 3, which has nontrivial solutions if and only if sin pa = 0, where a = d^ = a^ — a2- 
In other words, 



P 



Tin 
a 



1,2,. 



and the corresponding potentials are given by 

( A^^O 



. irnix — a2) . , , 

Ay =k2 sm —- /ki,Lu(y, z, t) 



Az =k-3, sin 



as - a2 

■nn(x — a2) 

as - a2 



h^,uj{y,z,t) 



c kj^ , Tmix — a2) 



as - a2 



h^,Lj{y,z,t) 



in the region a2 < x < a^- This is precisely the type 2 discrete modes discussed in [ig]. Their contribution to the 



:,TM. I 



Casimir energy is E^^^' . Notice that they are totally different from the transverse modes of type II. 
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In FIG. |3l we plot the Casimir forces as a function of mass when the background medium has constant refractive 
index ni, — 1 and ri}, = 2. In the graph, we choose a — ti = tr=10mn. From the graph, we can see that the TE 
contribution and the TM contribution to the Casimir force are approximately the same when the refractive index of 
the background medium is one. Since the first TM contribution to the Casimir force is equal to the TE contribution 
to the Casimir force in this case, we find that the second TM contribution to the Casimir force is negligibly small. 
However when the background medium has refractive index nt ^ 2, the graph shows that there is a significant 
difference between the TE and TM contributions. In fact, when m gets larger, the total Casimir force is dominated 
by the TM contribution. 

VII. THE CASIMIR EFFECT ON A PAIR OF INFINITELY PERMEABLE PLATES 

In this section, we consider the limits of the Casimir energy and Casimir force when the plates become infinitely 
permeable, i.e., ni^iir -^ co. In this limit. 



qT,l,qT,r ^ OO; QL.l = QL.r > \l ~ + ^1 '■= 90' 



Aj = Ai^-l, AJI = A:i^-1, Aj" = Ai"^^^:^ U, — r-:=A, 
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^Ib^lb ^bl^bl ^rb^rb ^br^br 



0. 



From these, it is immediate to show that the TE contribution to the Casimir energy for infinitely permeable plates 
is given by (PT|) . the same as in the case of perfectly conducting plates. It follows that the TE contribution to the 
Casimir force is given by (|22p . which is always attractive. 

For the TM contribution to the Casimir energy, one can show that 



OO /"OO 



^™^^J, J lnW(C,fci)fc^dfcide, (32) 



where 



^= (1 - A2e-2got,) (1 _ ^26-290*^) 

/ ASn _ p-2gotiUl _ p-2qot,.^ 

^ '\ (l-A2e-29ot<)(i_A2e-29ot.) 



(33) 



Therefore, the TM contribution to the Casimir energy Eq^^^ can be split into the sum of two terms E^'^l,' ^ and E^^^'' . 
The first TM contribution to the Casimir energy E^^^' is the same as the TE contribution, given by (PT|) . Therefore 
the corresponding Casimir force is also always attractive. The second TM contribution to the Casimir energy and 
Casimir force are given respectively by 

i^™^ " -A r r 1, (r A^(i-e-2.»*o(i-e-^^»--) ^_,.. A 

^^^ "4^27^ ]^ in(^l (^_^2e-2got,)(l„A2e-2'?o*.)^ jk^dk^d^, 

fTM, II ^ _ _^ / / J ji A e j ji A e j zg^.g 

^^ 2^2 y^ j^ 9L,b| ^2(i_g-2got,)(l_e-29ot.) 

When the refractive index of the background Uh is a positive constant, qL,b > Qo- Therefore < A < 1. It follows 

TM TT 

that the second TM contribution to the Casimir force F^^^ ' is also always attractive. This implies that the Casimir 
force acting on a pair of infinitely permeable plates is always attractive. 
Notice that 




90 



{ql.b - kl) - QL.b {ql - kl) ^ ^1^ + kle{iO - ^/m)^+^ 



90 



{qlb-kl)+qL,b{q^O-kl) ^^ + kliilO + ^m)^ + kl 
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where 



e(^) = 1 



m^c^ 



eb{uj)fii,ii^)h?uj^ 



Comparing to (|20|) . we find that the second TM contribution to the Casimir energy ([34|) can be identified with the 
TM contribution to the Casimir energy of a pair of dielectric plates with permittivity e that is due to a massless 
vector field. 
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FIG. 4: The dependence of the Casimir forces on the mass m when the background medium has refractive index 1 and 2. 

In FIG. m we plot the Casimir forces as a function of mass when the background medium has constant refractive 
index n\y — \ and riy, — 2. In this graph, we choose a = ti = t^^lOnm. We can see there is a significant difference 
between the TE and TM contributions to the Casimir force. 

VIII. THE CASIMIR EFFECT BETWEEN A PERFECTLY CONDUCTING PLATE AND AN 

INFINITELY PERMEABLE PLATE 



In this section, we consider the limits of the Casimir energy and Casimir force in Boyer's setup [23|, i.e., one plate 
is perfectly conducting and one plate is infinitely permeable. Without loss of generality, assume that the plate on the 
left is perfectly conducting, and the plate on the right is infinitely permeable, i.e., ei,/ir — >■ oo. In this limit. 
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The TE contribution to the Casimir energy and Casimir force is then given respectively by 

h 



pTE 



TE ^ "- 
Cas 27r2 
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/ / InM +e-29T,5(«,fci)aU^^^^^^^ 

qT.b{^,k±) 



oo poo 



exp(2gT,fc(^, fc_L)a) + l 



k±dk±d$,. 



(35) 
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It follows that the TE contribution to the Casimir force is always repulsive. 
For the TM contribution, one can show that 



^Cas — 



47r2 



OO /"OO 



1 w(e,fc^) , ,, ,, 



(36) 



where 



Wo 



=(1 - A^e-^**,) ([1 _ A]2 - [A + A]2e-29o*-) 
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(37) 



In this case, the TM contribution cannot be split into the sum of two contributions even if we assume that the 
background medium have unity refractive index. 
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FIG. 5: The dependence of the Casimir forces on the mass m when the background medium has refractive index 1 and 2. 

In FIG. m we plot the Casimir forces as a function of mass when the background medium has constant refractive 
index n;, = 1 and n\, = 2. In this graph, we choose a = ti = ir=10nm. From the graph, we find that the Casimir 
force is repulsive when the background medium has refractive index nb = 1. However, when the background medium 
has refractive index Ub = 2, the Casimir force can change from repulsive to attractive when the mass increases. This 
shows that the change of mass can change the sign of the Casimir force. 



IX. CONCLUSION 



In this article, we have derived the Casimir energy and Casimir force acting on two parallel plates due to the 
vacuum fluctuations of a massive vector field. We assume that the two parallel plates are made of real materials 
and they are placed in a magnetodielectric background. With Maxwell's equations replaced by Proca equations, we 
argue that the correct boundary conditions for massive vector fields are the continuities of (j), A, dxA^, Hy , Em , Bj^ and 



dtD 



2 2 



The last one is equivalent to the continuity of Dj^ in the massless case. Not all the boundary 



^ft2 ^ ^ 

conditions are independent. A set of independent boundary conditions is given by the continuities of A, H|| and the 
continuity of either </> or dxA^- The plane waves in an unbounded media can be divided into transverse waves of type I 
and type II and longitudinal waves, where the transverse waves of type I and type II are natural extensions of TE waves 
and TM waves in the massless case. For a system of several plane parallel layers of general magnetodielectric media, 
there are transverse modes of type I that satisfy all the boundary conditions. The contribution to the Casimir energy 
from these modes is the natural generalization of the TE contribution in the massless case, and we also call it TE 
contribution. In general, there are no transverse modes of type II or longitudinal modes that satisfy all the boundary 
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conditions. Therefore these two types of modes have to be combined. The formula of their contribution to the Casimir 
energy is very comphcated. However, in the massless hmit, the two polarizations reduce to one polarization and the 
TM contribution to the Casimir energy in the massless case is reproduced. Therefore, we also call the contribution 
from the superposition of type II transverse modes and longitudinal modes the TM contribution. 

In the limit where the plates become perfectly conducting, obtained by letting the permittivities of the plates 
tend to infinity, we write down explicitly the formulas for the TE and TM contributions to the Casimir energy. For 
general magnetodielectric background, the TM contribution cannot be split into two pieces. However, when the 
refractive index of the background is equal to unity, which happens for instance when the background is vacuum, 
the TM contribution can be split into two pieces, corresponding to the contribution from t ype 2 discrete modes and 
contribution from type 3 continuum modes discussed in the work of Barton and Dombey [l8|. In other words, our 
result reproduces the result of Barton and Dombey [31 for the special case they have considered. In this case, the 
Casimir force is always attractive. 

In the limit where the plates become infinitely permeable, obtained by letting the permeabilities of the plates tend 
to infinity, the TM contribution to the Casimir energy can always be split into a sum of two terms. In this case, the 
Casimir force is always attractive when the background has constant refractive index. 

To look for possible scenario that leads to repulsive force, we consider the configuration proposed by Boyer |23|, 
where one plate is perfectly conducting and one plate is infinitely permeable. In this case, we find that the TM 
contribution to the Casimir force cannot be split into a sum of two terms even when the background medium has 
unity refractive index. Numerical computation shows that the Casimir force can become attractive for nonzero masses 
when the refractive index of the background medium is not one. This shows that compare to the massless case, the 
Casimir effect of a massive vector field can behave very differently. 
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